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Abstract 

The case of a non-biased (light excited) p-n photo-voltaic (PV) cell working 
in a closed circuit with a resistive load is investigated by means of irreversible 
thermodynamics in the linear range. To our knowledge, this is the first time 
that both electronic processes, such as photo-excitations, carrier diffusions and 
recombinations, and electrical processes, such as current generation and the 
Joule effect, are considered together in the context of thermodynamics. An 
equation for the built-up voltage is obtained where the recombination resistance 
of the diode and a collection factor of the photo-carriers, depending on the 
thickness of the illuminated cell-side, are included in a natural way. 

PACS: 05.70.L; 72.40; 73.40; 85.30;84.60.j 



I. INTRODUCTION 

Photovoltaic (PV) cells are devices capable of photon-to-electrical energy 
conversion. To date, the most diffuse PV devices are based on (mainly silicon) 
p-n junction. The promise of new high efficient devices (dye sensitized, thermo- 
photovoltaic, hot carrier, organic) has reinforced the interest of researchers in 
this matter 1 . Besides the underlying chemical and physical processes and the 
technological features, since the seminal work by Shockley and Quiesser 2 , the 
main concern about the utilization of PV devices is solar energy conversion 
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and the related upper theoretical limit of conversion efficiency. The efforts 
made in succeeding in this task brought the attention of researchers in the 
thermodynamics of heat exchange between the black body source (sun) and the 
heat engine (converter). Most of the works that addressed this problem were 
grounded on reversible thermodynamics 3 . But improvements were obtained by 
considering the entropy production of PV conversion 4,5 . Recently, the efficiency 
problem was addressed by describing solar energy conversion as an explicitly 
irreversible process 6 . To this end, Onsager's reciprocity theorem was used to 
investigate near-equilibrium heat exchange. Despite the aforementioned efforts, 
literature lacks a study that explicitly takes into account the thermodynamics 
of electronic processes occurring in excited materials (carrier photo-generation, 
recombination and diffusion) as well as that of electrical current generation and 
the joule effect in a resistive load connecting the two PV-cell sides. For this 
reason, it is advisable to dwelling upon this topics since thermodynamics has 
an intrinsic generality so that in principle more differently chemical-physical 
based PV systems can be covered by the same algorithm implemented for a 
specific case. Differently from previously cited works, we are not interested in 
the efficiency problem. Thus, the paper is structured as follows. In Section 
2 some essential features of irreversible thermodynamics are briefly recalled as 
well as the main results of its application to the p-n diode. In Section 3, the 
theory is applied to PV cells and, thus, to PV closed circuits. The case of the 
biased PV cell is briefly considered. 



II. ENTROPY PRODUCTION IN A p-n DIODE 
A. General 

Before to beginning our investigation of the PV cell, it is best to clarify 
some points about irreversible thermodynamics. Let us first consider a (closed) 
system which can only exchange energy with its environment. Let this system 
also host chemical reactions that convert energies of exciting photons to heat, 
the state of reactions being characterized by the actual composition mixture 
of photo-sensitive molecules. Any change in the mixture composition causes 
a change of an associated entropy S c which will be referred to as " chemical" . 
Now, let some exciting photons cause departure from the equilibrium compo- 
sition. We suppose that the whole energy of absorbed photons is used only 
to change the composition mixture. This process is accompanied by a change 
5S C < 0. The excited mixture now returns to the ground state by reaching its 
original equilibrium composition after an entropy production SS' C > so that 
5S C + SS' C = . But the (dissipative) reactions bringing the mixture to the equi- 
librium composition also cause heat production. If the system were thermally 
isolated, the whole conversion process would cause its warming accompanied by 
an increase of its entropy. If the system temperature is maintained constant by 
heat exchange with a thermal reservoir, its entropy remains unchanged while 
increases the entropy of the thermal reservoir. In the stationary noncquilibrium 
state the system entropy (when it retains a meaning) is independent of time. 
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The case of an open system is more complicated. To explain this point, let us 
suppose that our system is now capable of releasing excited molecules to its 
surroundings and of receiving ground state molecules in such a way matter con- 
servation is allowed (we are interested in stationary states). Owing to escaping 
excited molecules, the negative change of the chemical entropy is not compen- 
sated by the positive one due to the dissipative reactions within the system. 
We also have that the heat produced by the dissipative reactions is released 
directly to the system's surroundings. Thus, apparently we have a net negative 
entropy production in the system. Nevertheless, the composition mixture is 
time-independent so that the chemical entropy does not change. In reality, the 
chemical entropy does not care about how the ground state is reached. Actu- 
ally, positive chemical entropy production is to be considered due to the flow of 
matter. In the range of linear phenomenological laws, the entropy production 
of chemical processes dS c /dt, has a minimum (compatibly with the system con- 
straints) in the stationary state (dS c /dt = 0) according to the minimum rate 
of entropy production (MREP) criterion 7 . Thus, if we are interested only in 
determining the current flow of the exchanged matter, the stationary condition 
could be sufficient. Instead, we could find a valuable tool in variational calculus 
if we are interested in characterizing the displacement from equilibrium of the 
composition mixture. 

To formalize the previous considerations, let us to write the configurational 
change of entropy as the internal entropy change diS as in textbook 7 



where is the chemical potential of the k-component of a mixture containing c 
chemical species and dn^ is the corresponding molar change. If at any instant 
the entropy changes as a function of chemical composition as well as other 
quantities characterizing the system, it is possible to write an equation for the 
rate of entropy production, that is, 



In eq. @ we take into account that dnk = v kid£,i , d£i stands for the 
displacement of the 1-th reaction, Vki for the stoichiometric coefficients , vi = 
d^i/dt, for the flux or velocity of the Z-reaction and T/ for the affinities, that 
is 7 , Ti = — ~^2 k VkiHk- It is assumed that near equilibrium the fluxes are 
linear with respect to the affinities, that is, vi = Li m T m , where L\ m are 
called the phenomenological coefficients 7 . As shown by Onsager 8 , based on the 
time reversal invariance of (microscopic) mechanical laws, the phenomenological 
coefficients form a symmetric matrix, that is Li m = L m i. In the case of photo- 
induced change of the composition mixture we have 




(1) 




(2) 



r 




(3) 
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where $ stands for the flux of the whole absorbed photons (in suitable units) 
and Cfc for the fraction of photon flux inducing molar change of the k-component. 
It follows from eqs that eq. © is to be replaced by d l S/dt = — $r$/T + 
J2i Tivi/T where T$ = J2k c k^k- 

Formally, the previous considerations about chemical reactions can be ex- 
tended to electronic processes in semiconductors: electron-hole pair photo- 
generations play the role of molecule photo-excitations; change of electron and 
hole level populations plays the role of composition change of chemical mixture; 
electron-hole recombinations play the role of dissipative reactions. Recombi- 
nations can generate photons and heat. We assume that the semiconductor 
temperature is independent of time and uniform. Of course, in the cases of PV 
cells inserted in closed circuits (open PV systems) part of the photon energies 
(that used in photo-excitation processes) can be lost in the connecting load as 
Joule heat or as some other electrical works. The flow of electron in the circuit 
allows to maintain unchanged the internal entropy in the stationary state. 

B. Entropy production in a p-n diode 

Near equilibrium (but also in most practical cases) the rate equations of elec- 
tronic processes can be written in the approximation of non-degenerate statis- 
tics. At the stationary state, densities of conduction band electrons and valence 
band holes depend on recombination processes and are independent of trapping 
in metastable levels 9 . Recombinations of conduction band electrons and valence 
band holes may occur directly, that is, by band-to-band processes, or indirectly, 
that is, by trapping in localized recombination levels. In any case, the departure 
of free carrier densities from equilibrium values depends on the affinity 10,11 

r = Hp + Mn = kT In ( = ] , (4) 
\npj 

n and p standing for the equilibrium densities. In the quasi-Fermi level (QFL) 
formalism, it is 11,12 ' 13 

r = F n - F p (5) 

where F n and F p stand for the QFL of free electrons and holes respectively. It 
can be shown that the rate of all the direct or indirect recombination processes 
can be written as v = LeqT, where Leq is a phenomenological coefficient 
accounting for all the processes involved, and that 11 Tivi = vT. Transport 
processes driven by electric and diffusion forces contribute to the entropy pro- 
duction depending on the current of electron ( j n ) and hole ( j p ) carriers, that 
is, 



diS Jn , j P , e j n ■ E e j p ■ E vT 
~ d T =a = — f" 9 rad M« - — ■ grad fi p h + — (6) 

where fj, n and [i v stand for the chemical potential of free electrons and holes and 
E = —grad V for the electric field. In eq. (jSJ we suppose that the temperature 
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is uniform in the semiconductor. The carrier currents are 

j n = —D n grad n - n x n E , (7) 

j p = —Dp grad p + p x p E . (8) 

where x n and x p stand for electron and hole mobilities respectively. At the 
stationary state, the continuity equations for the free carriers are 

div j n = -v , (9) 



div j p = —v . (10) 

that is, in the linear range, 

D n V 2 n-nx n V 2 V -v = (11) 

D p V 2 p + pXp\7 2 V -v = (12) 

It is to be pointed out that the continuity equations I jlll and I12|) have been 
proved to be consistent with the MREP criterion by considering only band-to- 
band recombinations 10 . The same calculations can be repeated by including 
all the recombination processes whose entropy production has the simple form 
vT/T as well as by including photo-excitations 



III. ENTROPY PRODUCTION IN A PV-CELL 
A. Differential equation for the affinity 

For simplicity's sake, in the following we refer to the diode scheme by 
Shocklcy 14 . Densities and field depend on the x coordinate only, the edge of 
the p-side is at x — and xt p and XTn mark the transition region in the p- 
and n-side, respectively. By using the Einstein relations, that is, eD p /x p — 
eD n / x n = kT, and by taking into account that in the approximation dealt 
with d 2 T/dx 2 w kT [(1/ra) d 2 n/dx 2 + (1/p) d 2 p/dx 2 ] , we obtain from the lin- 
ear combination of eqs l|llf> and (|12|) 

d 2 T , „ L ( D D D 



, kT^^ ( =2 + ni ) r = . (13) 
dx 2 D p D n \n p ' 



Let us define the free carriers lifetimes as 



P 



r n = , (14) 

kT Leq kT Leq 
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We introduce the characteristic length I of the affinity function defined by 



1 _ urr Leq ( Dp , D n 



- =kT ^_,^ + ^\ = - + _ 1 _ (15) 
V D p D n \ n p J DpTp D n T n 

so that eq. H13[) becomes 

d 2 T r 

As boundary condition for the solution of eq. I|ltj|) (in the case of forward bias) 
we require only that far from the junction the affinity vanishes (owing to the 
equality between electron and hole QFLs 14 ). Consequently, the affinity function 
is 



r (x) = r (x Tp ) exp [{x - x Tp ) /I] x < x tp , (17) 

r (x) = r (Xm) exp (x Tn - X/I) X > X tn ■ 

We suppose that the transition region is so small that xt p ~ xx n ~ Sp , where 
s p stands for the thickness of the p-region, so that T (xt p ) — T (xtu) = V a 
where V a is the voltage of the biased junction 14 . 

B. The case of photo-excitation with carrier diffusion 

In the case of band-to-band semiconductor excitation, any absorbed photons 
cause the generation of one electron-hole pair so that Ck = 1 (generalization to 
the multiple pair generation case is simple) and thus T$ = T. Equation 
including photo-excitation becomes 

In the approximation dealt with the continuity eqs. are 

D ^ n -^n^V-(v-^)=0 , (19) 

Equations (|19fl and l|20l) can be combined to give 

d 2 T 1 / $ 



dx 2 I 2 \ L 



J EQ 



. (21) 



Note that in the linear approximation, that is, by putting u = kT [p — p)/p~ 
r, eq. H21|l becomes similar to eq. (10) of ref. [15] (in which a general relation- 
ship between dark carrier distribution and photo-carrier collection in solar cell 
is presented) provided p is independent of spatial coordinates. Now, we suppose 
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that monochromatic radiation is used to excite the PV-cell at the edge of the 
p-side (x = 0). Thus, eq. (|2~T|) can be solved by assuming that the $ (x) is ruled 
by the Lambert-Beer law, that is, 

= ® exp(-x/d) , (22) 

d standing for the characteristic penetration length of radiations. Actually, 
d can be of the same order of I so that we should account for of the different 
properties of the solution of eq. (|21() in the cases of d smaller, larger or equal to I 
(this is important if one is interested in investigating the wavelenght dependence 
of PV efficiency). To simplify calculations, let us assume that d << I. With 
this approximation the solution of eqs. I|21|l can be written as 



T(x) = T + -nr- exp (—x/l)--pr- exp (-x/d)+V a exp (- \x - s p \ /I) 

\ l z L E qJ l l L EQ 

(23) 

where, in general, Tq and V a are not independent of each other (now V a cannot be 
considered the junction voltage). For x = we have T (0) = Tq + V a exp (—s p /l) 
. Of course, if s p » I only the interaction with photons allow for T (0) ^ 0. 
We point out that the value of Tq is fixed by the actual net rate of carrier 
generation in the PV-cell, that is, 

G = J {<P~v)d 3 x = AL EQ J - T^j dx 

A being the cross section area of the diode. In the case of the open PV-circuit, 
at the steady state G — so that (diS / dt) pv = holds. To better illustrate this 
point, we assume that the parameters of the affinity function are independent 
of x (in reality they can change when passing from the p- to n-region as well 
as near semiconductor surface). Moreover, we assume that the thickness of the 
diode (but not of the p-region) is much larger than I . Thus, by putting G; = Al 
and 0d = Ad we obtain (terms of the order {d/l) 2 are disregarded) 

G = (6 d $ - QiLeqTo) . (24) 

Now, it is easy to prove that at the open circuit condition j n — j p so that, by 
taking into account that (//„ + fi p ) = kT \i\np/rvp = T, 

¥)„~^+f-¥. 

whose integral over the whole PV cell volume is (see Appendix ) 

Zpv = [ {dS t /dt) pv d 3 x = ir G = . (26) 

Jcell 1 

C. Entropy production in the case of closed photo-voltaic circuit 
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Let an electric load be matched to the two diode sides, and the PV-cell 
operate in photo- voltaic mode (no applied bias). In this fraction of 

photo-generated holes sustains the electric current I ext so that the rate of re- 
combinations in the semiconductor turns out to be less than that of photo- 
generations (thus we are dealing with an open system). In this condition, the 
carriers diffusing towards the junction contain an excess of electrons with re- 
spect to the holes, such that it allows the generation of I ex t- On this ground, 
the external current is 

I ex t = G = e (e d $ - QiLeqT ) (27) 

The condition I ext ^ implies that negative and positive terms in eq. i|18|) are 
now imbalanced. To make this clear it is convenient to write the equation of 
entropy production in such a way that the external current will appear explicitly. 
Thus, by taking into account that 

3n-3v= I -fx ' (28) 

it follows from eq. I|18fl 

diS\ = _i_ . ar vT_^r_ li^d^ ( 

dt ) Pv T Jp dx T T T eA dx ' [ ' 

Form the integration of eq. I|29[) we obtain (see Appendix ) 

E PV = [ {dSi/dt) pv d 3 x + E En = -i— T(x Tn ) + E c (30) 

Jcell 1 e 

where we also considered the contribution 2c due to the carrier flow by the 
external circuit. Another way to obtain eq. (|30|l is to consider the entropy 
change due to current I ex t crossing the junction against the built up PV- 
cell voltage (charge separation in the photovoltaic process) V = T (iy n ) = 
F n (xtu) — Fp {xtu) 1A ■ The field at the junction is E — gradV so that (see eq. 
[T5)l . — \ ~j n — j pj ■ E JT = (I ext /Ae) dV/dx/T which after integration gives 
— (Iext/Te) V . Now, since T (xxn) ~ T (s p ) ~ To exp {—s p /l) it follows from eq. 

I ext T exp(-s p /l) 
-pv = -^j, 1" ■ (31) 

As for Sc we may calculate it by taking into account that electrons leave the 
n-side at the chemical potential (~p n ) n and holes leave the p-side (electrons 
from the conductor) at the chemical potential (~p p ) p so that, analogously to 
recombination processes, the associated affinity can be written as 



A 



(Mj„ + (-p P ) p (32) 
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which is nothing but the difference in QFLs of the two cell sides. Thus 



r-, A lext loo\ 

~c = f — (33) 

The affinity A is actually determined by the constraints due to the connected 
external circuit. In the case of a ourely resistive circuit A = e(R + R s ) I ext 
where R stands for the resistance of the connected load and R s for a resistance 
accounting for interface effects at the load contacts. 



r exp(-s p /Z) A I ext 

— JpTf + — 



"PV- 7fT^ iext + 7^, ■ (34) 

el 1 e 

Note that RI 2 xt /e also corresponds to the energy lost in the load (Joule heat) 
by the PV system . Now, by taking into account of eq. I|27(l we obtain at the 
stationary state 



j = eQd$o cxp(-s p /Q 

ext R + R s +exp(- Sp /l)/e 2 eiL E Q e 2 QiL EQ 

Equation lj35J) can be written in a simplest form by putting 



I g = eQ d $ , (36) 



Thus, we have 



Ro = J! ; W _ 7 = cxp (s p /l) g£ = cxp (- s p /l) R EQ , (37) 



Iext ' R + Rs + Ro l9 (38) 
As for the meaning of Rq it is convenient to rewrite it as 

Ro = exp (-s p /l) = exp (s p /l) R EQ 

Peq = 1/Leq being the recombination resistance (note that here the symbol is 
different with respect to the ones used for electrical resistances) 11,12 . R E q can 
be defined as the electrical recombination resistance normalized with respect 
to the volume having a thickness equal to the diffusion length. Now, the open 
circuit voltage turns out to be (R = oo) 

Voo = R I g = exp (-s p /0 R EQ I g (39) 

Actually, eq. I|37|l shows that if the thickness s p is very large with respect to the 
diffusion length, resistance Rq vanishes. This is a consequence of the increase 
of the time available for recombinations of photo-generated carriers before they 
reach the junction. 
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D. The case of a biased PV cell 



The case of biased junction is more concealed since T and V a are not inde- 
pendent. Indeed, by using eq. J22Jl, eq. (|27|) is to be replaced by 



Iext = e<dd$o-e@iL E QTo-2 [1 - exp {-s p /l)] e 2 <diL E QVa-e 2 <diL E QVa exp {-s p /l) 

(40) 

which only fixes a relation between the unknown variables I ext , Tq and V a . 
The stationary condition can only fix an additional relation thus leaving under- 
determined the problem. About the eq. Q4U|). note that as the thickness of the 
p-region becomes much larger than the diffusion length, all the photo-generated 
carriers recombine before the junction is reached. Thus, in this case, the elec- 
trical current is due only to the applied bias (diode mode), that is (note that 
the sign convention used for the current in the PV cells is the opposite with 
respect to the one used for diodes), I ex t = —2e 2 QiLEQV a — —2V a j Req — 
—I s (eVa/kT) which in the linear range is the current of the biased diode where 
I s = 2kT/eR,EQ is the saturation current 11 . We conjecture that, in the lin- 
ear range, V a = when Tq can match the external constraints, that is, when 
To exp (—s p /l) = \% + A where Vt> stands for the voltage of the biasing battery. 
Thus, at the condition I ext ~ we have Vb = RqI 9 = V^. 



Conclusions 

We investigated the PV conversion by means of irreversible thermodynamics. 
This has been made by considering explicitly photo-excitation, recombination 
and diffusion of carriers in the non-biased PV-cell as well as electric current and 
the Joule effect in a resistive load closing the PV-circuit. An equation for built 
up voltage and current in a closed PV-circuit is obtained where a collection 
factor and a recombination resistance are included. 



APPENDIX 

Equation l|26(l can be proved after direct and tedious calculations by taking 
into account that, in the case of non-biased PV cell (V a = 0), eq. (|19|) can be 
simplified as 

D v^P = (v-$) (41) 
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where Sp = p — p (note that, at the equilibrium, the change of carrier densities 
at the junction are connected to the built up field). Ea (|41l) which (see eq. 12311 
has the solution 

' - {(* + ££) ^ (' + £) ' } 

(42) 

In the case of open circuit, all photo-generated carriers recombine in the semi- 
conductor (G = 0) so that f cell {dSi/dt) pv d 3 x = 0. In the case of a closed 
circuit (Lext 7^ 0) we obtain from integration of eq. (|29|) 

J ^ {dS t /dt) PV d 3 x = ir (QiLeqTo - e d <t> ) - J ^d 3 x = 

where \i n (oo) is the chemical potential of electrons at the end-side of the n- 
region. By taking into account that in the p-region \x v « ~p p w /i p (s p ) (holes 
are majority carriers), \i n (s p ) w (oo) = — /x p (oo) (equality of QFLs at the 
edge of n-side 15 ) and that /i p (s p ) — \x v (oo) « r (s p ) (in the n-region affin- 
ity is determined by the diffusion process of holes through the junction), we 
have n n (oo) - jU„ (0) = {[/i„ (oo) + ^ p (s p )] - [/i„ (0) + [i v (s p )}} = T (s p ) - T . 
Therefore 

J (dS i /dt) PV d 3 x = ~ I -^T(s p ) . 

which represents only a part of the entropy production in the PV-cell since it 
lacks the contribution due to the carrier flow in the circuit. 
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